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ABSTRACT 

I constrain a possible supernova origin for gamma-ray bursts by modeling the dynam- 
ical interaction between a relativistic jet and a stellar envelope surrounding it. The 
delay in observer's time introduced by the jet traversing the envelope should not be 
long compared to the duration of gamma-ray emission; also, the jet should not be 
swallowed by a spherical explosion it powers. The only stellar progenitors that com- 
fortably satisfy these constraints, if one assumes that jets move ballistically within 
their host stars, are compact carbon-oxygen or helium post-Wolf-Rayet stars (type Ic 
or lb supernovae); type II supernovae are ruled out. Notably, very massive stars do 
not appear capable of producing the observed bursts at any redshift unless the stellar 
envelope is stripped prior to collapse. The presence of a dense stellar wind places an 
upper limit on the Lorentz factor of the jet in the internal shock model; however, this 
constraint may be evaded if the wind is swept forward by a photon precursor. Shock 
breakout and cocoon blowout are considered individually; neither presents a likely 
source of precursors for cosmological GRBs. 

These envelope constraints could conceivably be circumvented if jets are laterally 
pressure-confined while traversing the outer stellar envelope. If so, jets responsible 
for observed GRBs must either have been launched from a region several hundred 
kilometers wide, or have mixed with envelope material as they travel. A phase of 
pressure confinement and mixing would imprint correlations among jets that may 
explain observed GRB variability-luminosity and lag- luminosity correlations. 

Key words: gamma rays: bursts, supernovae: general, shock waves, relativity 



1 INTRODUCTION 

There is presently a growing body of circumstantial evi- 
dence linking some long-duration gamma ray bursts (GRBs) 
with afterglows to the explosions of massive stars. Super- 



in seven afterglows (iGalama et al. 1998; Bloom et alj 


199C 


Reichard Il998l: iTuratto et al.l I200Q: 


..Biornsson et al.l 


2001 


Lazzati et alJl200l|) although some 


of these could be light 



terglows of six other GRBs have been interpreted in terms of 
a wind-like ambient mediu m as expected around a massive 
star a t the end of its l i fe JPrail et alJll999l: iHoUand et al " 



2001 



iFrail et all I2OO0I: iLi fc ChevaheJ I2OOII: iPrice et al 



2002 h. However these are often equally well explained by a 
coUimated flow in a uniform medium. X-ray lines have been 
detected with moderate confidence in ab out half the after- 
glows for which they were investigated iToshio fc Daisuk3 
|2000|); however, the da ta analysis has been questioned (see 
iRutledee fc Sakdl2003l) . If real, these are most easily ex- 



plained by dense material su rrounding the burst engi ne, sug- 
gesting a stellar origin (e.g.. iBottcher fc FrveillioOlll . 

iFrail et'all i200lll have recently derived beaming an- 
gles for a number of GRBs from observations of their af- 
terglows. These authors derive gamma-ray energies rem- 
iniscent of supernovae: roughly 3 x 10^" erg for the 
observed lobes of GRBs. Similar results wer e reported 
by [ Panaitescu fc Kumar. (.200lh and iFreedman fc WaxmanI 
J2001D . 

Clear evidence that GRBs occur very close to mas- 
sive star formation would be almost as conclusive as a SN 
signature in an individual GRB. Several GRB afterglows 
show evidence for high c olumn densities (980703 and 980329; 
iGalama fc Wiierj|200ll) or high local gas densities (000926 
and 980519: iHarrison et alJ^OOlFWang et al.l2000l) . both of 
which connote star-forming regions. Likewise, the intrinsic 
extinction of GRB 000926 is characteristic of a gala x y disk 
JPrice et alj 1200 jl . | Bloom. Kulkarni. fc Diorgovskil j2003) 
have shown that the observed locations within hosts imply 
a tight correlation between GRBs and stellar populations. 
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considered too tight iBloom et al.lll999d) to be matched by 
merging neutron stars. Note however that the locations of 
merging neutron star pairs depends on their uncertain dis- 
tribution of natal kicks. 

If GRBs are a rare byproduct of star formation, rapidly 
star-forming galaxies should be over-represented as GRB 
hosts. In optical light host galaxies tend to look ordinary 
compared to contemporaries in the Hubble Deep Field, but 
[Ne III] and [O II] and infrar ed observations often in dicate el- 
evated star formation rates ijPiorgovski et all2001^ . At least 
eight afterglows h ave been associated with starburst or in- 
teracting galaxies ('Piorgovski e t alJl200ll : iFrail et al] 120021 : 
[cha rv. Bccklin. & Armus 2002). 

Although the association between long-duration GRBs 
and SNe is tentative (and applies only to the long-duration 
bursts for which afterglows are observed), the above evi- 
dence warrants a car eful evaluation. There are two ways a 
SN can create a GRB. Iciplgatel lll968h predicted that gamma 
rays might be produced in the very fastest, outermost ejecta 
of an ordinary su pernova explosion. This p ropo sal was re- 
cently revived bv iMatzner fc McKe3 j200Clll and iTan et alJ 
J200lf) . These authors showed that the GRB (980425) most 
compellingly associated with a SN (1998bw) is likely to be 
the result of trans- relativistic SN ejecta colliding w ith a stel- 
lar wind (see also iMatziier. Tan, fc McKed 1200 J) . In their 
model, as conjectured bv llwMncyto et al.l I^1998^ ~SN 1998bw 
was spherically symmetric or mildly asymmetric, and pro- 
duced the GRB in an external shock. Scaled-up versions 
of this model could produce external shock GRBs, at the 
expense of vast amounts (~ 10^* erg) of energy in nonrela- 
tivistic ejecta. 

In contrast. ISari fc PiranI (^§3) have argued that most 
GRBs require internal emission (e.g., by internal shocks) 
within unste ady ultrarelativistic win ds or jets, as originally 
suggested bv lR.ees fc Meszarod l)l994h . The arguments for in- 
ternal emission are strongest for rapidly-fluctuating cosmo- 
logical bursts wi th hard spectra, those least resembling GRB 
980425; also, see lPermer fc MitmanI l)l999h for arguments in 
support of external shocks. I shall assume for the purposes 
of this investigation that cosmological GRBs involve inter- 
nal emission within optically thin, ultrarelativistic outflows. 
For this to result from a SN, a jet must emanate from a 
star's core and pierce its envelope - shedding the baryons 
in its path - prior to producing the gamma rays observed 
at E arth. Such a ietlike expl osion is the conventional model 
fe.g.. lMacFadven fc Woosie'viil999i) for a supernova origin of 
cosmological GRBs. 

The goal of this paper will be to develop analytical mod- 
els for the phase of this latter model in which a jet, already 
created by the stellar core, must traverse the envelope and 
shove aside material in its path. These models, which ar e 
complementary to numerical simulations ijAlov et alJl200cll . 
are meant to elucidate under what conditions the hypothesis 
of a stellar origin is viable for the observed GRBs. In §f|21and 
|21l assume that jets travel ballistically within their stars; this 
allows one to place strict constraints on stellar envelopes. In 
^this assumption is reconsidered. It is shown that a phase 
in which the jet is hot, pressure conflned, and mixing with 
its environs would have interesting consequences. 



1.1 Stellar Progenitors 

Figure sketches the typical masses and radii of the stellar 
GRB progenitor candidates considered in this paper. In gen- 
eral, those that retain an outer envelope (e.g., supergiants) 
have quite large radii {R > 10 -R0) at the time of core 
collapse, whereas those depleted in hydrogen due to winds 
or binary interaction (e.g., those that have been through a 
Wolf-Rayet phase) are quite compact {R < 10 Rq). Among 
post- Wolf- Rayet stars, those containing helium ("He Wolf- 
Rayets" on the plot) are less compact than their He-depleted 
peers ("C/O Wolf-Rayets"). 

Very massive objects (VMOs) might have formed at 
high redshift due to the difliculty of cooling in the absence 
of metals and thus the large Jeans mass in primordial gas 
iHartouist fc Camero nll973). VMOs may also form today in 
rare conditions. Those initially more massive than ~ 25OM0 
die when their cores collap se to black holes, a nd are candi- 
dates for producing GRBs; iBond et alJ (^3) discuss their 
evolution. Again, the pre-coUapse radii of these stars depend 
on t heir mass loss. If p resent, their H envelope is quite diffuse 
(see iFrver et alJl2001^ . Otherwise the remaining convective 
core is very compact at the point of collapse. In addition 
to winds and binary mass transfer, VMOs may shed their 
envelopes in a super-Eddington phase during helium core 
contraction. Lone VMOs are thought to retain H envelopes 
if they are formed from sufficiently low-metallicity gas (al- 
though this is uncertain; Bon d et al. 1984:). 

The "supranova" model of lVietri fc Stellal posits 
that a SN explosion produces a rapidly spinning neutron 
star massive enough to collapse after sh edding its angular 
momentum jBaumgarte fc Shapirdll999ft along with ~ 10 
erg of rotational energy. The object collapses to a black hole 
with an accretion torus, firing a jet through the pulsar's wind 
nebula and its sheath of stellar ejecta. Because of the short 
viscous times of such tori, this model is most appropriate for 
short GRBs rather than the long bursts for which there is 
evidence of a SN connection. Also, we will see in ij^lthat the 
expanding stellar ejecta would prevent the production of a 
GRB entirely, unless it has either become Thompson opti- 
cally thin or has been cleared aside by the breakout of the 
pulsar wind nebula. The latter is possible in this model, be- 
cause the pulsar spindown energy exceeds the typical kinetic 
energy of SN ejecta. 

iFrver fc WooslevI Jl998h and lZhang fc Frveil (|20o3) dis- 
cuss a scenario in which the stellar core collapse leading 
to a GRB results from the coalescence of a h elium Wolf- 
Rayet star with a compact companion. Similar Iv. lFrver et alJ 
ll999 !) argue that binary mergers are likely to dominate the 
production stars whose cores collapse to black holes accret- 
ing through a disk. Progenitors created in this fashion will 
typically be stripped of their outer envelopes, hence com- 
pact. However the stripped envelope poses a potential bar- 
rier to jet propagation, as in the supranova model, unless 
it is strictly conflned to the equator of the system or has 
expanded to the point of being optically thin. 

In sections l2 . 1 l and fG . 1 I I examine the properties of stellar 
cores at the point of collapse. For this I shall assume that 
collapse sets in at t he oxygen ignitio n temperature, Tg = 
T/IO" K ~ 3.2 (e.g.,|Bondl^aL||l983). 
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Figure 1. Schematic of tlie typical masses and radii of the stellar progenitors considered in this paper f jLll . The plotted regions each 
represent several pre-supernova stellar models from the literature. 



2 JETS WITHIN STARS AND STELLAR 
WINDS 

Consider the progress of a relativistic jet outward from a 
star's core through the stellar envelope. Schematically, three 
distinct regions develop: the propagating jet; the head of the 
jet, where jet material impacts the stellar envelope, and a 
cocoon consisting of shocked jet and shocked ambient mate- 
rial. The se are familiar components from the theory of radio 
galaxies llBegelman fc CioffiiilQSS) : see figure |5| 

2.1 Presence of an uncollapsed envelope 

Stellar envelopes pose a problem for the propagation of GRB 
jets only if they have not collapsed prior to the launching of 
the jet. The stellar envelope's collapse timescale is greater 
(probably by a factor of at least a few) than its free-fall 
time iff. In figureOl compare stellar free-fall times with the 
intrinsic durations of GRBs for several possible redshifts. 
Only the very densest progenitors, the cores of very massive 
objects (iff = 23(Mcore/[100 M©])^/"* s at oxygen ignition) 
and helium depleted post-Wolf-Rayet stars (iff > 50 s) could 
plausibly collapse entirely in the durations of the longest 
GRBs, and then only at low redshift. In all other cases a 
stellar envelope remains to impede GRB jets. 

2.2 Head Motion 

For stellar core collapse to successfully produce a cosmolog- 
ical GRB, it must emit a jet that clears the stellar enve- 
lope from the observer's line of sight to the core. This is 



required if the jet is to achieve high Lorentz factors and if 
its internal shocks are to be unobscured by overlying mate- 
rial. This cannot occur if the Lorentz factor of the jet head, 
Ffi = (1-/3^)^^''^ (where fihC is the head's velocity), exceeds 
the inverse beaming angle XjQ. If it did, then the jet head 
would be causally disconnected from its edges and would be - 
have like a spherical blastwave jBlandford fc McKedllOTil . 
In this case essentially no material would escape sideways 
to form the cocoon. If instead Vh < i/6 then there is am- 
ple opportunity for shocked jet and envelope material to 
flow sideways and inflate the cocoon, so a successful GRB 
requires 

Fh < 1/6 (stellar envelope). (1) 

Observers of GRB afterglows often identify an achromatic 
break in the light curve with the d eceleration of the swept 
up sheU from T^y 1/9 to Th <l/9 (iRhoadJigQST) . In order 
for this to be possible, the head must have satisfled 

Ffi > 1/S (circumstellar medium) (2) 

while it was still being driven forward by the jet (prior to 
its deceleration) . Condition ||5J only applies to GRBs whose 
afterglows exhibit such a break in their afterglows. 

The expansion velocity of the jet head is given to a very 
good approximation by the balance of jet and ambient ram 
pressures (momentum fluxes) in the frame of the jet head. 
This approximation is most accurate when ambient material 
is cast aside into the cocoon, as is the case if equation 10 is 
satisfled. Ram pressure balance means 

Pjhj(TI3)]hC +Pj = Paha{TI3)lc^ + Pa (3) 
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Figure 3. GRB durations compared to stellar dynamical times. The histogram of observed burst durations (t9o/(l + z), for redshifts 
2 = 0, 2 = 1 and 2 = 5) from the current BATSE catalog I http://cossc.gsfc.nasa.gov/batse/ I is plotted to show the timescales for long 
and short bursts. Notable bursts with afterglow observations are plotted as diamonds and labeled. Also plotted are the typical free-fall 
times for carbon-oxygen ("WR: C/O") and helium ("WR: He") progenitors evolved f rom Wolf-Rayets, blue and red supergiants ("BSG" 
and "RSG," respectively), the primordial 300 Mq progenitor ("VMO") calculated bv lFrver et al J 1200 J) , and the IOOMq core of a VMO 
with no envelope. For GRBs arising from SNe, this plot shows that an uncoUapsed stellar envelope remains to be traversed by the jet — 
except in the case of especially long GRBs coming from compact carbon-oxygen stars or VMO cores. 
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where {ri3)jh is the relative four- velocity between the jet 
and its head, {Tf5)h is the four- velocity of the head into the 
ambient medium, p, p, and h = {e ~\-p)/{pc?) are density, 
pressure, and enthalpy (e is total comoving energy density), 
and the subscripts j and a refer to jet and ambient material. 
As Pa ^ paC? in a stellar envelope /la = 1 to a good approx- 
imation and Pa can be ignored on the right-hand side. To 
leading order in T"^, pj can be ignored on the left. For a 
stationary ambient medium, these approximations give 

Pi 



jMarti et alJll994l) where 



We! 

dM/dr ' 



(4) 



(5) 



The second equality, which holds to leading order in T"^, 
is derived by noting that the kinetic plus internal energy 
density of the jet, evaluated in the lab frame, is T^hjpj — 
TjPj —Pi\ its isotropic luminosity Liso is this quantity times 
47rr^/3j ; and the mass per unit length of the ambient medium 
is dM/dr = Anr^pa- Here r is the radius and the true jet 
luminosity is L = (^^/4)Liso. 

The relation between Fh and Tj implied by equation 
||1J| takes simple limits in two regimes. If L > rt, then the 
rever se shock into the jet is nonrelativistic (|Sari fc PiranI 
I1995I) and Fj — Fh <^ Fj. This case is physically unattainable 
within a star and is achieved outside only if the ambient 
density is extremely low. In the opposite limit of a relativistic 
reverse shock. 



= 1/2 



2F|, Ft>L>l; 



I3h \ Ph, L « 1 ^ ' 

(see also iMeszaros fc Wa:x:ma'nll200ir i. Therefore, L is useful 
in determining the observer's time tobs/(l + z) = t — r/c for 
the jet head to expand to radius r if viewed at redshift z. 
Since dt — dr/{l3c), 

dtobs _ dr 
1 + z ^ UF'c 



(7) 



for both non-relativistic and relativistic jet heads, so long as 
L <C F*. In presupernova envelopes and stellar winds dM/dr 
is relatively constant, and since Liso is likely to vary slowly L 
can be approximated with its average value. More generally, 
one might know how Liso varies as a function of t (which 
is also tobs/(l + z) at the origin, r = 0). If one also knows 
M(r) (e.g, from a model of the star and its collapse), then 
equation l|7| integrates to 



_d£_ 

l + z 



I dr , 



/ dM 



(8) 



giving tobs{r) implicitly. 



2.3 Cocoon Structure 

The jet cocoon is the region containing spent jet material 
and shocked ambient material. Its extent is equal to that of 
the jet, but its width is determined either by pressure bal- 
ance with the surrounding gas, if there is time for this to be 
achieved, or else by the expansion of a lateral shock into the 
envelope. The latter case holds so long as it predicts a cocoon 



pressure Pc in excess of the hydrostatic pressure or collapse 
ram pressure pa in its environment. Let us first consider the 
case of an adiabatic cocoon whose pressure exceeds that of 
its surroundings, and check this assumption in H6.ll Let us 
also restrict attention to the case where the head velocity 
is subrelativistic (L < 1, eq. in which case the cocoon 
pressure Pc is roughly constant away from the jet head. 

The cocoon created by a jet of constant opening an- 
gle expands self-similarly so long as its width Rc expands in 
proportion to its extent r, and so long as no other size scales 
affect its structure. Under these conditions, numerical sim- 
ulations can determine coc oon structures exactly. H owever, 
analytical estimates (e.g.. iBegelman fc Cio ffi' 1989V while 
not as accurate, elucidate how cocoon properties scale with 
L and 9. 

The cocoon expands nonrelativistically in a direction 
normal to its surface at the velocity /3c given by 



PaC^f^c = Pc 



(9) 



If pa is evaluated at the point where the cocoon is widest, 
the normal direction is sideways and thus f3c — Rc/{ct). 
The pressure pc is related to the energy Ein deposited in the 
cocoon and the cocoon volume Vc through pc — Ein/ (SpaVc) 
(an overestimate, as part of _Ein is kinetic). If a cocoon of 
length r and width Rc is idealized as a cone, 14 — {-k /?i)R^r 
and pa — Pa(r) in equation @. But, since r = /3hC and 
Rc ~ /3cc, Vc ~ {Tr/3y{l3c/PH f. 

Now, Ein is the energy emitted in time to catch up with 
the jet head at radius r. (After breakout, Ein is available to 
drive an observable outflow and pos sibly a precursor: see 
iRamirez-Ruiz. Celotti. fc Reed i2002ll and Jll) As long as 
L ^ Fj the flight time of the jet can be neglected compared 
to tobs(7')/(l + z), so -Ein(r) is essentially all of the energy 
emitted up to then: 

/•tobsM/{l + z) 

Ein{r) = / L dt 

Jo 

4 1 1/2 



16 



LisorM{r)c 



(10) 



using equation (|7|l. Along with the expressions for pc, Vc and 
/3c, this implies 

^3/8^1/2^ 



(11) 



This implies that the cocoon is nonrelativistic (and slower 
than the jet head) so long as the jet head is also nonrel- 
ativistic, for L < 1 in this case (eq. |S]), and < 1 for a 
collimated jet. 

2.3.1 Cocoons vs. Jet-Driven Explosions 

Equation IIH assumes that the length of the cocoon is set 
by the advance of the jet head; this requires /3c < Ph. If 
this condition is violated, the cocoon expands around the 
jet and develops into a roughly spherical blastwave (see 
alsolgeEelman & Ciofh 1989). For a nonrelativistic jet head, 
equations © and (llll l indicate that /3c < so long as 
L>e*. For a more precise criterion, consider the velocity 
/3bw of a spherical blastwave powered by two jets: 



/3b 



0.90 



2Ein{t = r/Pb^ 
M{r)c^ 



1/2 



0.746' 



,2/3 f 1/3 



(12) 
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w here the coefficie nt is given by the PGA/K approximation 
of lOstriker fc McK cc ( 1988), for a wind bubble with a ratio 
of specific heats 7 = 4/3 in a medium with pa oc r~^. The 
head outruns this blastwave {f3h > /3bw) if L > (0/90°)'', 
i.e., 

%^>(<^)'*/Wc'. (13) 

This condition can only be violated if L < [9/90°)'^ < 1, 
so that only jets driving nonrelativistic heads can violate it. 
I assumed a nonrelativistic head in deriving equation 11311 : 
however there is nothing to suggest that relativistic jet heads 
can be swallowed by their cocoons. 

A successful GRB requires that stellar envelope mate- 
rial be cleared from the path of the jet so it can sustain 
Fj > 10^ at r ~ 1 AU where internal shocks are thought to 
form. A jet-cocoon structure is required for this: a spherical 
blastwave would not accomplish it. For this reason GRBs 
must satisfy (0/90°)'' < L < within the stellar envelope 
(eqs. 0, and 1131 ) so that a cocoon forms but does not 
overcome the jet. Of these limits on L, only the lower need 
be considered because L <^ 9~'^ for any reasonable combi- 
nation of opening angle and stellar model. 

One can use equation (1131 in equation (IIUI to eliminate 
either Liso^/c, which gives a lower limit on Ei-n{r), or to 
eliminate 9*M{r)c^, which gives an upper limit on Ein{r). 
In terms of the total energy per lobe Ein{R) deposited in 
the stellar envelope, the formation of a cocoon rather than 
a spherical blastwave implies 

O-lOe^MonvC^ < Ein{R) < 0.61Lisoi?/c. (14) 

A spherical blastwave is not sensitive to the site of energy 
injection, whether at a jet head or near the collapsing core, 
and for this reason the upper bound applies to energy from 
any source that is entrained in the stellar envelope before the 
jet breaks out of the star. Note that the lower bound implies 
that Ein{R) must be at least a fraction (0/37°)^ of the rest 
energy of the envelope in the jet's path, (0^/4)MenvC^. 



2.4 Timing Constraint 

One further and very important constraint derives from the 
durations of GRBs. In the internal-shock model, fluctuating 
gamma-ray emission reflects variability in the central engine 
and p ersists only while this engine is running jSari fc PiranI 
Il995f) . Define t^^^bs as the observed duration of a GRB, in- 
cluding any precursor but excluding its afterglow. In the 
coUapsar model, this must be preceded by the cocoon phase 
during which the jet crosses the stellar envelope. The central 
engine must therefore be active for at least ^,obs + tohsiR) 
in the observer's frame. It is unlikely - though not impos- 
sible - for tj^obs to be much shorter than tohsiR), ^o^' this 
would require the central engine to shut off just as its effects 
become observable. So 

^4^<W, e,E,4R)<E„ (15) 

where e-y is the efficiency with which the jet's kinetic lumi- 
nosity is converted into gamma rays in the observed band. 
The second inequality derives from the first as long as the 
jet's luminosity is relatively constant. 



2.5 Jet variability 

I have assumed a constant luminosity jet up to this point, 
but GRBs are often observed to fluctuate significantly in 
intensity on very short timescales. How should the above re- 
sults be adjusted for jet variability? First, note that the pro- 
cess producing gamma rays (e.g., internal shocks) is likely 
to accentuate the intrinsic variability of the source. Second, 
the observed propagation speed of the jet head, dr/dtoba, 
is equal to cL^^^ (eq. |7|). A fiuctuating jet thus progresses 
more slowly than a steady jet of the same mean luminosity. 
If one uses the average value of L to constrain a star's mass 
and radius by requiring that tohs{R) is not long compared 
to the observed burst (eqs. 1151 above and 1181 below) , then 
this constraint is only tightened if one accounts for variabil- 
ity. Similarly, a variable jet deposits more energy per radius 
than a steady jet of the same mean luminosity; this only 
tightens the constraints derived by requiring a jet-cocoon 
structure (eq. 



3 CONSTRAINTS ON STELLAR HOSTS 

To apply the above constraints to stellar progenitors for 
GRBs, observed quantities must be related to the the pa- 
rameters of equations 1131 . 1141 . and (1151 . This is possible 
for cases where, in addition to tj^oha, the redshift z and jet 
opening angle 6 have been derived from afterglow observa- 
tions. To construct L one requires the isotropic (equivalent) 
kinetic luminosity Liso and the mass per unit radius dM/dr 
in the environment. 

I adopt for GRBs' energy, luminosity, and duration 
the following definitions: L-yiso = ^?7,iso(l + z)/tj^obs', and 
{E-y, L-y} = {E-y^iBo, ^7180 } X 6^/4:. HeTS tfie isotropic gamma 
ray energy is -E^.iao, and is the average efficiency factor 
relating the gamma-ray energy to the total (kinetic, Poynt- 
ing, and photon) luminosity of the jet. The net energy E 
and luminosity L represent only the approaching jet, which 
presumably has a counterjet. 

A key assumption employed throughout this section is 
that the 7-ray half opening angle 6 is equal to the jet's half 
opening angle while it crosses the outer stellar envelope. 
This is essentially the same as the assumption that the jet 
is ballistic (rather than hot and pressure-confined) at that 
point; this assumption is revisited in ^ I assume the v alue 
of 9 derived from afterglow observations jRhoadjliggll can 
be used to characterize the jet as it crosses the outer stellar 
envelope and emerges from the star. The above definitions 
identify L^iso as the mean value (energy E^ in source-frame 
duration t7,obs/(l + z)). This is somewhat arbitrary, since 
GRBs are highly variable; see H2.5l for justification. 

For numerical evaluations I use either L^iso.si = 
L^iao/lO^^ erg/s, or £7,50.5 = E-^/lO^"-^ erg. The former is 
a characteristic value for GRBs and can be obs erved with- 
out d etermining 6. The latter is motivated by iFrail etlil] 
J200ll) 's result that E^ = IQ^O-^^O-^ erg for ten GRBs whose 
6 could be determined. 

In a stellar envelope the radial average value of dM/dr is 
Mbuv/R, the ratio of envelope mass to stellar radius. Indeed, 
many pres upernova envelo pes have density profiles close to 
Pa oc r~^ ^Chevaliejll989^ ■ for which dM/dr is constant at 
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its average value. The average value of L is 

L = 1.30x 10-^^«iAJ^ (envelope). (16) 

In order for the jet head to be relativistic (L > 1), the star 
must have a mass per length much lower than Mq/Rq, or 
the gamma-ray efficiency Ej must be small. 



have esti mated d by r equiri ng that — IQ^"'^ erg as sug- 
gested bv lFrail et alJ ll200j) . These points fall in a narrow 
band on the plot because Ljieo/O* = (1 -I- z)Ej/{6'^tobs)- 
With E-y held fixed, the dispersion in this quantity is domi- 
nated by 6^'^ (2 dex rms, for the points plotted) rather than 
(1 + z)/tohs (0.5 dex rms). 



3.1 Constraints from Burst Duration and Cocoon 
Formation 

Under the assumption of a ballistic jet, the observed time of 
jet breakout is 

U^^iR)=MA{l + z)(-^^^\^\. (17) 



7iso,51 



This, along with the constraint tobs(i?) ^ t7,obs (eq. 1151 ') il- 
lustrates that typical long-d uration GRBs are most easily 
produced in compact stars l|MacFadven et alJ 120011^ . This 
constraint is best expressed 



Menv R ^ ^ 



^7,obB/(l -I- Z) 

64.4 s 



(18) 



the left-hand side pertains to a hypothetical stellar model 
and to the efficiency of gamma radiation and is constrained 
by observables on the right. Note that the combination 
i7iao,5i [t-,,oha/{i + z)]'^ is related to an observed burst's flu- 
ence and duration through its comoving distance (i.e., phys- 
ical distance at redshift zero), rather than its luminosity 
distance. Note also that the above constraint is independent 
of 0. 

In figure 0] I illustrate the derivation of I/isotobs/ (1 + ■^)^ 
and constraint 11811 from the observed fluence and dura- 
tion and an observed or estimated redshift (for Qm =0.3, 
A = 0.7, Ho = 65 km s~^Mpc~^). Bursts with well-observed 
redshifts are plotted as solid d iamonds; those for which 
iFenimore fc Ramirez-Rui3 ll200Cf) have estimated redshifts 
from a luminosity- variability correlation are plotted as open 
circles. A progenitor model passes the timing constraint if 
it lies to the left of its burst. Most of the bursts plotted 
are consistent with e-yMBuvR'^WMQRQ; for none of them 
does this limit exceed QOMqRq. This constraint is nearly 
independent of redshift and derives primarily from the dis- 
tribution of (fluence • tobs) among GRBs. VMO cores and 
Wolf-Rayet stars are compatible with many bursts. Blue su- 
pergiants are compatible only with the very brightest and 
longest. Red supergiants and VMOs that retain their en- 
velopes are ruled out, as are optically-thick shrouds of ex- 
panding ejecta which might persist in the "supranova" or 
stellar-merger models of central engines. 

Equation (I13II gives a complimentary constraint on the 
basis that a jet-cocoon structure exists: 



Mcnv Rq j 



17.0° 



(19) 



This constraint is derived from observations of bursts in fig- 
ure |S| Similar to figure |1] filled diamonds here represent 
bursts for which afterglow observations have all owed an es- 
timat e of 6 in addition to Liso (as reported bv lFrail et alJ 
I2OOII . and ref erences therein): open circles repre sent redshifts 
estimated bv IFenimore fc Ramirez-Ruid ll2000f) . For these I 



4 CONSTRAINTS ON STELLAR WINDS 

GRBs from supernovae are likely to occur within a dense 
stellar wind, an environment t hat potentially affects both 
the GRB itself and its afterglo w ( Kumar fc Panaitescul2o'ool : 
iRamirez-Ruiz. Drav. Madau. fc Tout 20011 . 

The internal-shock model for GRB emission posits that 
significant variability on a time scale 5t (= 5tobs/(l + z)) 
arises from the collision of shells within the jet at radii of 
about 

ns = 2V]c St, (20) 

as discussed by iRees fc Meszaro'3 il994l) . ISari fc PirarJ 

fl995) argue that an external shock with the ambient 
medium - to which the afterglow is attributed - cannot 
create a fluctuating gamma ray burst. For internal shocks 
to occur, the external shock (jet head) must move beyond 
rIS on a timescale not long compared to the duration of the 
bur st: AtobB(ns) <A.obs . 

ISari fc PiranI l|l99!fl define Np = t7,obs/<5fobs as the 
number of pulses that fit within the burst duration given 
a characteristic separation Stohs- With equations JSJ, 12UII . 
and this definition, the criterion Atobs(ns) J; i7,obs becomes 



NS ViooJ I iV„ 



(21) 



(22) 



laari fc PiranI argue that Np ~ 100 is typical, altho ugh a wide 
varie ty of time scales is observed within bursts ijLee et alJ 
I2OOO!). 

For a wind dM/dr is the ratio of mass loss rate to wind 
velocity, M^/v w 5 SO 



^ 5.89 X 10" L^iso-oiVw.s , . 
L = — : (wmd) 



(23) 



where Vw,i = 
M„/(10-5 Me yr 
ues for Wolf-Rayet stars 



/(lO* cm s"^) and M^u-s = 
are normalized to c haract eristic val- 

B), 



l|Nugis fc LameriJl200m . although 
the pre supernova m ass loss rates of such stars is not well 
known iLangedll99ll) . Expressing condition (|22|l in terms of 



(24) 



r,<35of^V'' I^Wi^^^/^ 

VlOOy V e-,M^,-s 

The analogous criterion for a uniform ambient medium was 
presented by^iran ( 1999). This upper limit on Fj must ex- 
ceed the lower limit required for the escap e of the observed 
gamma rays, described most recently by iLithwick fc Saril 
(2001). These authors inferred Fj > 340 for five of the ten 
ordinary bursts in their Table 2; equation 1241 suggests that 
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Figure 4. The product of a GRB's fluence and duration constrains the product of mass and radius in any stellar envelope through which 
the GRB jet (if moving ballistically: see JHJ could have emerged rapidly compared to the period of emission. The observed quantity 
fluence-toi,s a^nd the redshift imply the quantity Lisot^j^_,/(l + z)^ intrinsic to the burst. By eq. 1181 . this places an upper limit on the 
model parameter e-^Ms^^R, where e-y is the efficiency of gamma-ray production; so, a model should lie to the left of its burst on the 
plot. This limit insensitiv e to GRB redshift when ^ > 1. In this plot solid diamonds are the bursts plotted in figure with redshifts 
and fluences derived from lpYail et all i200 J) and references therein; circles are 220 bursts for which iFenimore &: Ramirez- Rui3 i200Cl) 
estimate redshifts using a luminosity- variability correlation (these authors did not allow z > 11). The burst illustrated with dotted lines 
is compatible with a carbon-oxygen post-Wolf-Rayet progenitor if e-, < 15%, or with a helium-bearing post-Wolf-Rayet star or the 100 
Mq bare core of a very massive object (VMO) if e-y ^3%. (For other VMO core masses, use MenvR oc M^/"^ at O ignition.) It is not 
compatible with blue or red supergiant pr ogenitors (BSGs or RS Gs), VMOs that have not lost their radiative envelopes, or the expanding 
envelope of the "supranova" progenitor of lVietrF fc Stellal il999D (not shown), which moves to the right on this plot. Similar conclusions 
can be drawn for almost all of the bursts plotted. Lines of constant z arc shown for A = 0.7, Q,m = 0.3, Hq = 65 km s~^ Mpc^^. 



these could not have occurred in a stellar wind representa- 
tive of Wolf-Rayet stars - unless, for instance, Np ^ 100 
(which does not appear to characterize the bursts in their 
Table 2). 

However, condition 11241 might be circumvented by the 
efficient sweeping-forward of the amb ient medium by run- 
away pair production, as discussed bv lMadau fc Thompson' 
1 20001) ■ Thompson & Madau (200(1) and [Bcloborodov 
I 2002h . If the relativistic flow has generated a fraction of 
the observed photons at a radius smaller than ris, then these 
can clear optically-thin ambient gas from the region. 

A spreading break (albeit one difficult to detect: 



iKumar fc Panaitescul l2000l) will occur in the afterglow if 
equation ^ is satisfied. The remaining isotropic kinetic lu- 
minosity after the gamma- rays have been emitted is (1 — 
e-y)L^iBo/£T Along with 1231 . this becomes 



> 1.64° 



(1 — e-^)L~fiso,5iv-w,s 



1/4 



(25) 



This criterion could only be violated by the narrowest of 
GRB jets in especially dense stellar winds. 
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Figure 5. The opening angle and isotropic luminosity of a GRB - in the combination Liso9~*, eq. I19i — constrain the ratio of mass 
and radius in any envelope through which the GRB's jet can travel without being engulfed within its own bubble. (The jet is assumed to 
be moving ballistically; see As in figure |1] a viable model should lie to the left of its GR B on the plot. Her e, the constraint is most 
sensitive to 9. Filled diamonds repre sent bursts with inferred from afterglow observations iFrail et al.ll200ll and references therein). 
Open circles are the bursts for which IPenimore fc Ramirez- Rui3 i200Cl) estimate redshifts; for these d is estimated under the assumption 
of a common energy — 10 ' erg (causing 9 to correlate tightly with L^\^a). Most bursts are compatible with all stellar progenitors 
by this criterion; however, those with large 9 are not. Diffuse stars most easily pass this test. VMO cores fail first, when 9 ~ 10° (plotted 
is a 100 Mq core; M^nv/R oc M^/^ at O ignition). The burst illustrated by dotted lines is compatible with RSGs, VMOs that retain 
their envelopes, BSGs, and helium post-Wolf-Rayet stars for values of e-y less than 90%, 50%, 3%, and 1%, respectively. The supranova 
model travels to the left on this plot as it expands, and easily passes this constraint. 



5 BREAKOUT AND BLOWOUT: 

OBSERVATIONAL CONSEQUENCES 

The energy deposited by the jet as it crosses a surround- 
ing stellar envelope {Ein in is available to power (con- 
cievably) observable phenomena distinct from the GRB and 
its afterglow. These include transients associated with the 
emergence of the jet head from the stellar surface (break- 
out), and the escape of hot cocoon material along the path 
of the jet (blowout). Note that a type Ic hght curve is not a 
product of jet energy, as SN Ic are powered by radioactive 
nickel, which is not produced at the jet head. 

The energy Ein is typically less than the total GRB jet 
energy (for successful GRBs) according to condition 11511 . 
Nevertheless the effects of deposited energy may be of ob- 
servational interest if they cover a wider opening angle, dom- 
inate a different observed band, are emitted with higher ef- 
ficiency, or occur where the jet stalls before breakout. 



front followed by additional postshock acceleration iSakurail 
Il960l) . Sufficiently energetic explosions can produce rela- 
tivistic ejecta this way, so long as the stellar progenitor 
is sufficiently compact jMatzner fc McKedHooil) . The im- 
pact of these ejecta with a circumstellar wind may pro- 
duce a transient of hard photons; this is a plausible ex- 
planation for the association between supernova 1998bw 
and GRB 980425 jM atzncr fc McKedll999l:lTan et alj|200ll : 
iMatzner. Tan, fc McKeai200lh . ' 

iMacFadven fc WooslevI Jl999t). iMacFadven et all 
teOOl'). and lZhang. Wooslev. fc MacFadvenI ll2003l) have sug- 
gested that shock breako ut might lead to the X-ray precur- 
sors seen in some burst s l|Lar os et al ]|l984HMurakami et alJ 
| l99lt lin 't Zand et alJ Il999^ ■ IMacFadven fc WooslevI and 
iMacF^^er^^-l.! concentrate on the prompt flash from 
the shock-heated stellar photosphere; however, the fast 
ejecta have a gre a ter st ore of kinetic energy to be tapped. 
iLvutikov fc Usovl (l2000l) discuss several cases of precursor 



5.1 Shock Acceleration and Breakout 

In ordinary supernovae the fastest ejecta are produced at 
the surface of the star. There, the rapid decrease in stel- 
lar density leads to a whip-like acceleration of the shock 



1 SN 2003dh, which appeared in the afterglow of GRB 030329, 
appears similar to SN 1998bw; however, the Tan et al. model 
cannot explain a two-peak GRB of ~ 10^^'^ erg emerging from a 
SN of comparable energy. 
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activity in gamma rays, including one burst for which iron- 
line e mission has been suggested (GRB 991216: Ivietri et alJ 
l200lfl . If precursors are due to shock breakout then their 
energy is a small fraction of Ein{R) and is therefore 
limited by conditions 114t an d (115t; however it may yet be 
observable, lin 't Zand et al.l l)l999|) argue that the X-ray 
precursor of GRB 980519 resembles its X-ray afterglow to 
the extent that the afterglow effectively preceded the GRB; 
as argued by Pacz vnski ( 1998il . this is consistent with the 
emission from material with 7 < Fj ejected prior to the 
GRB emission. It is therefore worthwhile to estimate the 
distribution of ejecta kinetic energies from shock breakout 
in a jetlike explosion. 

To do this one must identify at what point the jet head 
makes a transition from the state of ram pressure balance 
described in H2.2I to the state of whip-like shock accelera- 
tion discussed by Tan ct al. (2001). The jet head obeys ram 
pressure balance only if shocked ambient material exits the 
jet head more rapidly than the head accelerates. While this 
holds, the forward ambient shock cannot travel far ahead 
of the jet reverse shock. At some point near the surface, 
however, ambient material cannot exit the jet head prior 
to breakout. The forward shock will then accelerate away 
down the density gradient, and t he flow becomes progres- 
sively more normal to the surface (|Tan et al.ll200ir i. 

If f}± is a typical value for the perpendicular component 
of velocity in the jet head, and x = (R—r) /Ris the fractional 
depth within the stellar envelope, then material is trapped 
within the head if {9R)/I3i_ > {xR)/f3h, i.e., if the time to 
exit the head exceeds the time to reach the surface. A rea- 
sonable guess for /3x is that shocked ambient gas exits the jet 
head at the postshock sound speed, in the frame of the head. 
For a nonrelativistic head, then, f3± ~ (8/49)^/^/3^. In the 
relativistic case, the transverse velocity saturates at ~ c/-\/3 
in the head's frame; in the star's frame, p± ~ l/{3^^^Th) 
jRhoadjll999l) . So the transition occurs when 



(49/8)1/2^, 
3^/"Fh<^. 



L{x) < 1; 
L{x) > 1 



(26) 



where {4>nr,4'r) ~ 1 are uncertain parameters. In the rela- 
tivistic case, x is determined implictly once Fhix) is known. 
For this, apply equation Q to the outer density distribution 
of the stellar progenitor 



Pa = Ph 



( — 



(27) 



where n is the effective polytropic index and the coefficient 
Ph is an extrapolation to r = R/2. In a radiative outer 
layer p^ can be derived fro m the mass, radius, and lum i- 
nosity of the pro genitor star (jKippenhahn fc Weiger1^[l994^ : 
iTan et al] ll200l[) present formulae (their equations [25], [45], 
and [48]) for p^ in Kramers or Thomson atmospheres. 

Once X and Phix) or Vhix) have been identified, t he 
production of fast ejecta follows from iTan et all i200ll) 's 
theory (their eq. [54]) if one matches the shock four- velocity 
7s/3s with the head four-velocity ThPh at the depth x. With 
equations I^J and 1251 1 setting a reference depth, external 
mass, and shock velocity, their theory predicts an isotropic- 
equivalent ejecta kinetic energy 



-1.65ri7pp-(l-5S7p-l) 



-Efc,iso(> F/) 



Lis 



2.687p- 



L\so 



R 



(28) 



when the transition occurs in the nonrelativistic regime, and 



-Efe,iso(> F/) 



{<t>re)~ 



OJBlT-rp 



^(1.587p-l) 



(n + 2.70) exp 



L'lSO 



R^'phc^ 



(n+3.26)(n+6.69) 
1.22n(n+4) 

4.327p -3 
7i(4-,p-3) D 

c 



(29) 



when the head is relativistic at the transition. Here, 7p = 1 + 
1/n is the envelope's polytropic exponent. Thesen equations 
have been simplified by the restriction Tf ^ 1. 

In general, the appropriate value of Ek.iso{> F/) is the 
minimum of the values given by equations 1281 and 12911 . 
Note that i5fc,iso(> F/) is an isotropic equivalent; the total 
kinetic energy per lobe in ejecta above F/ is smaller by 6^ /i. 
If one varies 9 holding the other quantities fixed, this total 
energy above F/ is maximized when the two expressions are 
equal, whereas the isotropic value continues to rise slowly as 
6 is decreased in the relativistic regime. In the relativistic 
regime, the above formulae only appy to F/ > F(,(a;)^'^^ - 
i.e., only to those ejecta involved in quasi-spherical shock 
acceleration. 

The uncertainty of {(pnr, 4>r) x 9 leads to a much greater 
uncertainty in Ek,iso{> F/) for the nonrelativistic than for 
the relativistic case: for instance, when n = 3, Ek,iso(> F/) 
varies as {4>nrd)~''''^ in the nonrelativistic and as {4>r9)~'^'^^ 
in the relativistic regime. Physically, this difference arises 
because ram pressure and shock acceleration give very dif- 
ferent velocity laws in the nonrelativistic regime {(3h oc Pa^^^ 



and Ps oc pa 



respectively) whereas they give very sim- 



ilar velocity laws in the relativistic regime {Th oc pa^^^ and 
Fs oc pa"'^"^^, respectively). Another uncertainty concerns 
whether the head velocity should be matched to the shock 
or postshock velocity; the latter choice increases the ejecta 
energies by only 6 — 7% for both relativistic and nonrela- 
tivistic transitions. 



5.1.1 An Example: the SN 1998bw Progenitor 

To illustrate these estimates of jet-driven shock breakout, let 
us consider the proge nitor model for SN 1998bw adopt ed by 
IWooslev et alj lll999t) and studied bv lTan et all j200lD. Th e 
parameters n, /?, and ph can be derived from iTan et all 's 
tables 2 and 3. Adopting their fit for the outermost re- 
gions (whose rest energy is < 1.7 x 10^'' erg), n = 4, 
ph = 335 g cm~^, and J? = 1.4 x 10^° cm. The isotropic 
kinetic energy in breakout ejecta is therefore: 



-Efe,iso(> F/) 



1.9 



7iso,51 



f<t>r9V 

V2.5° J 



2.1 



^7iso,51 



/ 4>nrO \ ■ 
V 2.5° / 



X 10^«F7«-^^5 erg. 



(30) 



(n -I- 2.70) exp (1.98 + 7.03/71 - 1.52n) 



If this formula were to give a result > 2 x 10^" erg, a different 
envelope fit would be appropriate; however, the qualitative 
result would be unchanged. 

Even in isotropic equivalent, the energy predicted by 
equation 13011 is small compared to a cosmological GRB. It is 
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comparable in magnitude to the values derived bv lTan et alJ 
for the spherical explosion of SN 1998bw. While sufficient to 
produce GRB 980425 at a redshift of 0.0085, it would not 
contribute to the appearance of a GRB at 2 ^ 1- 

It should not be surprising that the energy of motion in 
shock breakout is intrinsically much smaller than that avail- 
able in the jet, as the accelerating shock is powered by the 
jet for a brief period (a small range of radii) prior to break- 
out. Breakout does produce a spray of ejecta with a variety 
of Lorentz factors, which may produce weak transients if 
observed off the jet axis. 

5.2 Cocoon Blowout 

The cocoon inflated by a jet prior to breakout is filled with 
hot gas that is free to expand away fr om the star after break- 
out, constituting a "dirty" fireball (|Paczvnskilll998h which 
may be visible either through its thermal emission, through 
its circumstellar interaction, or through line absorption and 
fiuorescence. The cocoon energy Bin (eq.EU comprises 



4.4 X 10* 



-J7iso,51 



Mo, 



Mc- 



R 

Rq 



1/2 



— erg. 

3°/ ^ 



(31) 



iRamirez-R.uiz. Gelotti. fc Reej (|2Qq3) have considered 
observable implications of cocoon blowout under the hypoth- 
esis that no mixing occurs between shocked jet and shocked 
envelope material, so that the cocoon fireball can acceler- 
ate to Lorentz factors comparable to Tj . A likely alternative 
is some mixing between jet and envelope, especially near 
the jet head where shear is strong. Suppose the shocked jet 
mixes with the envelope material within an angle (f)mi^d of 
the jet axis before entering the cocoon: ^mix — 1 for mix- 
ing near the jet head only. The energy per mass of the hot 
cocoon material is then {2/ (j)rai^)L}^^ c? , corresponding to a 
subrelativistic outfiow velocity 



\/2 ?! 



0.27 

0mix 



/4 



Rg; I \M,nJ 



7130, 51 



1/4 



(32) 



where I assume L ^ 1 as is appropriate for the compact 
progenitors favored in [0 If blowout occurs into an angle 
Sout , then the expanding cocoon material becomes optically 
thin (r < 1) along its axis after roughly 



14 



^7iso,51 



(1) 



1/4 



0.4 cm^ g 



1/2 



{9/3°) 



3/4 



hours. 



(33) 



Before this occurs, however, photons will diffuse sideways 
out of the blowout cloud when the lateral optical depth 
riat ^ Ooutr and lateral velocity vut satisisfy riat ^ c/ (Svut)- 
This occurs at ~ 13 hours for the fiducial parameters listed 
above (independently of Sout); the amount of internal en- 
ergy persisting in the ejecta is quite small by that point, 
indicating that a thermal pulse will not be observed unless 
the fireball is enriched with short-lived isotopes. The modest 
cocoon energy is available for circumstellar interaction but 
may be masked by the GRB afterglow. The cocoon blowout 



does, however, represent a screen between the jet-afterglow 
shock on the outside and the expanding SN ejecta and cen- 
tral engine on the inside. 



6 PRESSURE CONFINEMENT 

Two assumptions made in sections |5| and |3 remain to be 
checked, both concerning the effects of an external pressure 
on the cocoon or on the jet. H6. II addresses the assumption 
that the cocoon drives a strong shock into the stellar enve- 
lope. H6.2l concerns the possibility that the jet is self-confined 
by its own cocoon pressure, then spreads sideways outside 
the star before producing gamma rays. 



6.1 Ambient pressure confinement 

If the ambient pressure is greater than the pressure within 
the cocoon, then the assumption of a strong shock in equa- 
tion @ is incorrect and the cocoon does not inflate as de- 
scribed in H2.3I The cocoon is overpressured relative to the 
envelope by (jScC)^ /c^, where Cs is the envelope's isothermal 
sound speed: 



c.(r)^ = ^ = a(r)-^^-^W 



3r 



(34) 



With a{r) so defined, the virial theorem stipulates (a) — 1 
when the mean is weighted by binding energy and (l/a) = 1 
when weighted by thermal energy within the star. In general, 
Q(r) ~ 1 wherever the scale height is of order the radius. (In 
a polytrope of index n, a is related to IChand rasckhar 1939|'s 
variable v by 3/q = {n + l)v.) Using equations IIH and I34II . 
Pc > p when 



361 



2/3 



\c) G 

GM{r) 



4a 



17.3 



7/3 ^5 



(35) 



This constraint is expressed in terms of the total jet luminos- 
ity L; I shall evaluate i t at the fiducial co llapse temperature 
Tg ~ 3.2 (O ignition: ' Fond et alj[l98i) . The appropriate 
value of 9j, though different from O-y, cannot exceed unity. 

In ordinary core-collapse supernovae (below the pair in- 
stability limit) the core is degenerate but the envelope above 
the collapsing core is not; also, gas pressure dominates over 
radiation pressure. At the oxygen ignition temperature, con- 
dition 1351 1 becomes 



L>3.2 X 10™a" 



-1 7/3 



erg s 



(36) 



where ^ is the mean mass per particle in a.m.u. 

In contrast to ordinary supernova cores, VMO cores are 
dominated by radiat ion and are appro ximately n = 3 poly- 
tropes in structure iBond et alJll983l . They obey i?corc = 
0.26[Afcorc/(100 Mq)]^/'^{3.2/Tq) Rq; hence pc > Pa when 



> 



1.3 X lO"a^/^0f ^ 
7/6 



Mc. 



100 Mq 



erg s 



(37) 
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The critical luminosities in equations 1361 and 
13711 should be compared to the value L ~ 3 x 
10'^°(1 0%/e-^)(10 s/t^) erg s~^ implied by iFrail et alJ 
||2001j)'s standard value E-y = 10^°-^ erg. In both ordinary 
supernovas and collapsing VMOs, the jet cocoon could pos- 
sibly be pressure-confined in the collapsing region (Tg > 3) , 
but not in the hydrostatic region where Tg < 1. 



the primary assumption being that the jet head expands 
nonrelativistically {L < 1). Setting pj = pc and using Tj = 
I/O-,, 

2'^{cB?j'f'-Kpa = O^r'^L (cocoon — confined, hot jet) (39) 

which provides a means to evaluate pressure confinement if 
one has a constraint on Rj. 



6.2 Self-confinement by cocoon pressure 

The constraints on stellar envelopes derived above assume 
that observational determinations of the intensity (Liso) and 
opening angle of gamma-ray emission can be applied to the 
earlier phase in which the GRB crossed the stellar envelope. 
This requires that the jet is dynamically cold in the stellar 
interior, so that it does not widen once outside the star. It 
is important to consider the alternative: that the jet may 
contain internal energy sufficient to spread in angle after 
breakout. In this case the emission angle 6~, will exceed the 
jet angles 9j attained within the stellar envelope. The jet 
would then be more intense within the star than outside, 
vitiating (or at least relaxing) the constraints on stellar en- 
velopes derived in §S|3and|3| 

For the flow to expand after breakout it must possess 
relativistic internal energy, hjpj ~ 4pj, and it must be slow 
enough to spread, Vj < l/6j. (Note that cold jets will be 
heated and decelerate locally to Tj < 1/9 j, if they are 
crossed by an stationary oblique shock.) In this case the 
final opening angle will be set by the relativistic beaming 
of the jet at breakout: 6^ ~ l/Fj; a suggestion made by R. 
Blandford. In general, 

9 J ^ fT-7- :r+dArconi) 

-L j t,^conf ) 

conf ) ■) (^conf )1 (38) 

where the two possibilities refer to hot, confined jets (with 
the confining pressure is released at rconf) and cold, ballistic 
jets, respectively. ^ One must have rconf — R for pressure 
confinement to be significant; if rconf ^ R, the ballistic-jet 
constraints of the previous section would still hold for the 
outer envelope. For confined jets, Fj (R) is determined by an 
observational determination of the final opening angle. At 
the same time, when hjPj ~ 4pj , and so long as Fj is still 
significantly above unity, the jet luminosity is 

L = 47r_RjFjPjC (hot jet) 

(see the discussion after equation [S|), where Rj = 9jr is the 
jet radius. If evaluated at rconf, Fj may be replaced with 
1/9-,. 

What determines Pj? The confining pressure is clearly 
time-dependent in this scenario, which complicates matters, 
but one can set pj ~ pc at rconf. The cocoon pressure is 
paiPccf' where /3c is given by equation (lllll : 

r3 3 9 r>2 4 4 

L pa = TV C RjT Pc, 

^ IZhang. Wooslev. &: MacFadverl l2003l) suggest eq. I38i does not 
describe their numerical results; however, it is not clear what 
rconf or Fj (rconf) is applicable to the wide-angle ejecta they see. 
Although the choice of a single rconf is perhaps oversimplified, it 
remains a useful parameterization. 



6.2.1 Adiabatic jets 

One such constraint arises if the jet expands adiabatically 
from its launching region, since this provides a formula for 
Rj (pj). For the steady flow of a relativistically hot gas {pj 3> 
pjC^) at relativistic speeds {(3j ~ 1) through a channel, the 
simultaneous conservation of luminosity L = (4F|pj)7r_RjC 
and mass flux AI = FjPjTrR^c, along with the relation pj oc 
pf^, imply 

rj(r) = — and Pj(r) ^ - — 5-F"'' (adiabatic jet) (40) 
ro 47rrgC ■' 

where ro is the radius from which the flow accelerates 
through Fj = 1 and should reflect the dimensions of the 
central engine, tens or hundreds of kilometers. In equation 
1391 1. this implies 




The 7-ray opening angle produced by pressure confinement 
of an adiabatic jet is thus lower than any of the derived 
opening angles, the discrepancy being much greater in terms 
of pressure (since Pj/pc oc 9~^'^ at fixed ro). 

A jet expanding adiabatically from ro ~ 10 km can 
therefore be confined in the core of the star, but would be- 
come cold and ballistic in the outer envelope. We found 
in JB] that the ambient core pressure may be high enough 
crush the cocoon, suggesting that the jet can also be con- 
fined by ambient pressu re in this reg ion. This scenario was 
described by^eszaros fc ReesI (120011) : as they note, it holds 
only in the stellar core and fails (producing a ballistic jet of 
fixed 9j =9-,) in the envelope. However, we must question 
whether ro might also change. 

6.2.2 Nonadiabatic jets: crossing shocks and mixing 

The above argument shows that a jet expanding from a 
launching region ro of order 10 km will no longer be pressure- 
confined in the outer stellar envelope. But, what if the prop- 
agation is not adiabatic? This can result from shocks that 
cross the jet, or from mixing between the jet and its environ- 
ment (cocoon or stellar envelope). Heuristically one expects 
the jet's memory of ro to be erased in either process; a larger 
effective ro would ameliorate the discrepancy highlighted in 
perhaps making pressure confinement more realistic. 
First, consider a stationary oblique shock that crosses 
the jet. The postshock jet pressure (pj2) is brought into 
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equilibrium with pc, the external pressure that launches 
the shock. The shock front must be causal with respect 
to the postshock flow, i.e., sound travels across the (post- 
shock) jet faster than the shock itself does. If the postshock 
Lorentz factor is Vj2, sound fronts can propagate at an angle 
l/{\/^Vj) to the jet axis. This must exceed the angle 9ca of 
the crossing shock relative to the jet axis (although usually 
not by much : Iwhit hami 1 9 74ll . In order for the shock to cross 
the jet, one must have O^s > 6j\ so, Y'j2 ^ l/(\/3^cs) < . 
The jet therefore makes a transition to a hot, pressure- 
confined state at its current radius; the new eflective value 
of ro for the postshock jet is Rj/Tj2 > \/i9csRj- 

One can, in fact, estimate and the shock obliquity 
9cs from the expression Lj = 4iTR^r^2Pj2C (valid for hot 
jets), from the condition < l/{V3dcs), and from the re- 
quirement that pj2 — pc- If Pc = plpaC? is estimated via 
equation mn . one finds 



3.5 



Menv Rq lO^^erg/ s ( 3 
Mq R Lj 



1/8 



(42) 



which can be read as an estimate of the number of cross- 
ing shocks along the jet axis. Note that eq. 1381 predicts 
9-y ~ 6cs, roughly speaking, if the jet is effectively released 
after one of its crossing shoc ks. These results appear broadly 
compa tible with the work of Zh ang. Wooslev. fc MacFadverJ 
i2003r . and indicate the potential importance of this route 
to pressure confinement. 

Second, consider the mixing of ambient material into 
the jet during a pressure-confined phase. Nonadiabatic mix- 
ing can be described by imagining relativistically hot jet 
material flowing though a uniform channel of fixed cross- 
sectional area Aj . At some initial time the channel is allowed 
to mix with some area dAj of ambient material (density pc, 
pressure pc) at rest. Assuming perfect mixing and nonrela- 
tivistic surroundings (pc <K pc(?), the change in jet param- 
eters Pj, pj, and Tj are determined by the conservation of 
energy, momentum, and rest mass through the mixing event. 
That is, 



^^jhjPj-Pj 
^]Pjh,p, 
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dA-i 



(43) 



Algebra then yields differential equations for pj , pj , and Tj 
as functions of A. 

A couple interesting conclusions can be drawn from this 
exercise. First, the heuristic argument that ro would be for- 
gotten is correct: the effective value of ro = Rj/Tj increases 
monotonically in non-adiabatic jet expansion, because Rj 
increases whereas Vj decreases (or if Pc ~ Pc = 0, stays 
constant). This in turn implies that pressure confinement is 
extended somewhat by mixing. 

Secondly, mixing forces the nondimensional jet param- 
eters to trace a characteristic trajectory that whose impli- 
cations would be observable. Recall that the final opening 
angle 6'-y is the inverse of Fj if the jet is confined. Also, final 
Lorentz factor F/ is equal to the ratio of the jet's luminosity 
to its mass fiux, 



1 F, 



In terms of these dimensionless parameters, equation 14311 



can be restated 
dF, 



dlnMj 
din M, 



= — F" 



(4F| 
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1) 



r7,(2F^ + l)+F3 



%+2F2(,,, 



2) + F,3 



%(2F +l)+F^-F, 



(44) 
(45) 



where dlnMj = pcdA/{T jpjA) represents the fractional 
mass per unit length mixed in. The ratio of these equa- 
tions (in which the denominators cancel) describes set mix- 
ing trajectories for rjjiVj), as shown in fig. 1^1 There exists 
an attractor solution, which is 77 = 4F| (i.e., pj — TjpjC?) 
when Fj ^ 1 and approaches (??j,Fj) (1, !)• Mixing tra- 
jectories rapidly join onto it by decreasing either Vj or T]j as 
mass is mixed in. 

As shown in figure |S| mixing trajectories do cross 
through reasonable values of F/ and 9^. Does this mean that 
these jet parameters arise from mixing? This seems plausible 

if there ever exists phase of pressure-confined evolution in 

1/2 

the stellar core or mantle. If 2Fj < 77 , adiabatic expansion 



and mixing both act to bring 2Fj 



1/2 



If 2Fj > r] 



1/2 



the 



adiabatic tendency for Fj to increase is counterbalanced by 
the nonadiabatic decrease in Fj. This makes the mixing at- 
tractor solution a probable one for the end state of jets that 
are pressure-confined at some stage. (A more detailed anal- 
ysis would require a dynamical analysis of jet expansion.) 

If so, then the attractor solution would imprint the cor- 
relation 



11° 



100 



1/2 



(mixing attractor) 



(46) 



which may have observational implications in that L^-^so oc 
V f for a fixed Lj/e^. Such a trend has been implicated in the 
correlations between spectral lag and luminosity, variability 
and luminosity, and afterglow break time an d luminosity 
among GRBs (e.g.. iRamirez-Ruiz fc Llovd-Ronninel l2002l : 
ISalmonson fc Gala!mall2002r) . The mixing attractor provides 
a motivation for a rather tight correlation among bursts with 
similar jet luminosities. 



6.2.3 Confined jets: timing properties 

The variability-luminosity correlation would be affected also 
by the filtering of jet variations by pressure gradients, which 
occurs only when the jet is pressure-confined. Now, the in- 
teral dynamical time of a relativistic jet is ~ Rj/c, or F j Rj / c 
in the lab frame. If a strobe approached the observer at 
Fj, pulsating once per dynamical time, the observed period 
would be foreshortened to 



Sto 



Ri 



c 



(47) 



The same result can be derived by considering a standing- 
wave pattern of wavelength Rj in the jet's frame, or Rj /Vj 
in the lab frame, that is swept past a fixed pressure-release 
radius rconf and leads to jet pulsations as the peaks and 
troughs pass by. As we saw above, crossing shocks and mix- 
ing increase the effective value of ro and thus increase the 
characteristic variability timescale. Shocks and mixing also 
decrease L^iso by decreasing Fj and thereby increasing 9^; 
this implies a correlation between jet variability and 7-ray 
brightness. 
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Figure 6. Evolution of pressure-confined jets due to mixing witli their surroundings (eqs. 1441 and 1451 '). The ratio r/j of luminosity to 
mass flux is plotted versus jet Lorentz factor; these set the final Lorentz factor Ff and opening angle 8-y. Breaks in the trajectories are 
spaced by a doubling of jet mass per unit length. Mixing causes jets to move down and to the left along the plotted trajectories toward 



Vj 



1. In contrast, adiabatic evolution causes jets to move horizontally rightward to higher Fj (smaller at fixed rjj and Ff. No 



solutions exist below the dotted line (rij = Tj 
rectangle. Inside this region, mixing correlates observables: higher Ff corresponds to lower 



and are estimated by the dotted 



6.3 Pressure Confinement and Envelope 
Constraints 

What should we make of the possibility that pressure- 
confined jets will evade the constraints on stellar envelopes 
derived in §321 and |3| for ballistic jets? Equation 14H shows 
that jet confinement works if 

a value that does not appear to violate timing constraints 
(eq. \47\ ). This indicates pressure-confined jets are plausible 
if shocks or jet mixing increase ro sufficiently. 

Why would a jet change from pressure-confined to bal- 
listic? This could occur in a couple ways. First, the jet could 
accelerate in such a way that sound no longer crosses it de- 
spite it being hot, i.e., rj(r) > 1/6 j(r) for some range of r. 
However this appears to require the ambient density to drop 
to a near vacuum, implying it would not occur spontaneously 
within a star. Alternatively, the jet could become internally 
nonrelativistic fe^PjC^)- As we have seen, adiabatic evo- 
lution works in this direction whereas crossing shocks and 
non-adiabatic mixing oppose it. The transition would then 
depend on the specifics of jet propagation and mixing. Fi- 
nally, the jet could be confined by internal magnetic stresses, 
a possibility highlighted by the 7-ray polarization discovered 
bv lCoburn fc Boggj (|200a), but beyond the scope of this pa- 
per. 



Jet pressure confinement would naturally lead to time- 
dependent jet properties, since the confining pressure evolves 
during the cocoon phase and after breakout. Crossing shocks 
and mixing with ambient material are also unlikely to be 
steady processes. One would expect evolution during GRBs 
as a result; however, no tre nds are discerni ble except for 
changes in pulse asymmetry jLee et alJl200ol) . 

Apart from this contraindication, it is difficult to judge 
whether pressure confinement persists across entire stellar 
radii. If so, then the constraints on stellar envelopes in §^11 
and 121 are relaxed. Numerical simulations can potentially 
solve detailed questions such as this one. However, caution 
should be used in interpreting them, because of the many 
decades separating the launching radius from the stellar ra- 
dius. 



7 CONCLUSIONS 

This paper constrains possible stellar progenitors for GRBs 
by requiring that jets of GRB-like luminosities and dura- 
tions can clear a path for themselves in the star's envelope 
prior to producing gamma rays. Of all the possible progeni- 
tors, only the compact carbon-oxygen post-Wolf-Rayet stars 
(SN Ic progenitors) and the bare cores of very massive ob- 
jects can plausibly collapse in the durations of long GRBs at 
low redshift ( t|2.1l figures and |3J. In more extended stars, 
the outer stellar envelope remains to impede the progress of 
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the jet, and it is likely to violate observational constraint. 
Therefore, I conclude that GRBs with SN progentiors come 
primarily from type Ic or VMO-core events. Type lb is not 
thoroughly ruled out by this work, but type II (supergiant 
stars) are. 

The most stringent constraint on a stellar envelope 
arises from the requirement that the jet can traverse it in 
an observed time not much longer than the duration of the 
GRB. Under the assumption of a ballistic rather than a 
pressure-confined jet, this constraint is independent of the 
inferred opening angle of the burst, and (given an observed 
fluence and duration) depends on its inferred comoving dis- 
tance rather than its luminosity distance. This makes it in- 
sensitive to an uncertainty in redshift. Given the luminosi- 
ties and durations of GRBs (regardless of their redshift; fig. 
[IJ, only post-Wolf- Rayet stars and VMO cores are compact 
enough to satisfy this cr iterion. The variability-lum i nosity 
corr elation discussed bv iFenimore fc Ramirez-Ruid (|200(J) 
and iReichart et al.l (1200111 allows this constraint to be ap- 
plied to a large number of bursts in the BATSE catalog 
(figure 0J|. Only a few are compatible with blue super- 
giant progenitors; red supergiants and VMOs with envelopes 
are ruled out. Also r uled out is the "supranova" model of 
IVietri fc Stelial lll999D (a supernova followed by a GRB) , un- 
less the SN ejecta are optically thin by the time of the GRB, 
or unless the pulsar nebula is energetic enough to clear them 
aside. 

Post-Wolf-Rayet stars have bee n favored among stellar 
GRB progenitors since the work of iMacFadven fc WooslevI 
il999t) . on the basis that their compact envelopes delay the 
GRB jet breakout the least. The above constraint quanti- 
fies and strengthens this conclusion, and relates it to the 
observed properties of GRBs rather than those of a specific 
model for the central engine. 

For ballistic jets, an additional constraint arises from 
the requirement that the GRB jet cocoon should not over- 
take its driving jet and produce a spherical explosion. This 
is generally not as restrictive as the constraint from burst 
durations discussed above, but it does become important 
for jets with opening angles exceeding 10° (figure I^J. VMO 
cores, being the most compact, are the progenitors most 
likely to be ruled out this way if they have not collapsed 
prior to the GRB. The requirement of a jet-cocoon struc- 
ture also gives interesting upper and lower bounds on the 
energy entrained in the stellar envelope during the phase 
of jet propagation (eq. 1141 '). This energy is store d in the jet 
cocoo n and is available to drive a "dirty" fireball jPaczvnskil 
^9^ of expanding cocoon material after the jet breaks out. 

The breakout phenomenon is itself a candidate for pro- 
ducing a transient, as may have happened in SN 1998bw to 
produce GRB 980425. In section|S]it is calculated how much 
kinetic energy is channeled into relativistic envelope ejecta 
during a jet's breakout, by matching the propagation law for 
the jet's terminal shock onto the relativist ic shock and post - 
shock acceleration behavior described bv lTan et alj i200ll) . 
I find the ene rgy of this ejecta to be small compared t o that 
of the burst. iRamirez-Ruiz. MacFadven. fc Lazzati ll2002ll 
have recently discussed how the upscattering of photons 
from the shocked envelope by the jet may produce a precur- 
sor of hard gamma rays; however, note that this is reduced 
in importance for the com pact stellar progenitors favored 
by the timing constraints. IZhang. Wooslev. fc MacFadvenI 



i2003h have suggested breakout ejecta as the origin of short, 
hard bursts; however their estimate of the energy is at odds 
with that calculated here. 

If the gamma-ray photons are not able to clear 
away a stellar wind in th e region around the star i n 
the manner describ ed by iMadau fc Tho mpsonI i200Ct) . 
iThomDson fc Mada^ J2000f) . and lBeloborodovl J200j) . then 
the presence of this wind places an upper limit on the jet 
Lorentz factor. This limit arises in the internal shock model 
for GRB emission because the presence of the external shock 
limits the distance within which internal shocks can form. 
The equivalent limit has been pr esented previously for uni- 
form ambient media lPiranI 1^93); however, for a sufficiently 
dense stellar wind it can conflict with the l ower limits on jet 
Lorentz factor (e.g.. iLithwick fc Sarill200lh . 

Figure |7| illustrates the above criteria for the specific 
case of GRB 000418, assuming a ballistic jet with = 10%. 
Its seven-second duration is briefer even than the free-fall 
times of VMO cores. B ecause of its large opening angle 
(11°; iBerger et all 1200 ll) . it would not succesfuUy form a 
jet-cocoon structure in any uncoUapsed portion of the VMO 
core. In fact, a jet of its inferred luminosity could cross noth- 
ing more extended than the most compact of Wolf-Rayet 
stars in the GRB duration. I conclude from this that it came 
either from a compact carbon-oxygen Wolf-Rayet star, or 
from a VMO core that managed to produce a GRB of briefer 
duration than its free-fall time, or that it did not have a su- 
pernova origin. These restrictions change quantitatively, but 
not qualitatively, if e-y ^ 10%. 

Pressure confinement of jets provides a way for GRB 
hosts to evade the constraints on stellar envelopes listed 
above for ballistic jets. This occurs because, for a given 7-ray 
opening angle, confined jets are narrower and more intense 
than their ballistic counterparts. The computational results 
of Alov ot al. ( 200j) and [Zhang. Wooslev. fc MacFadve^ 
(2003) correspond to pressure-confined jets. In Hfci.2.1l it was 
shown that jets that do not mix w ith their environs canno t 
remain pressure-confined (see also iMeszaros fc ReedlioOlT) . 
Simulations that do not resolve the jets' launching scale may 
not observe this effect, however. Mixing of jet and envelope 
{ W).2.2l is capable of sustaining a pressure-confined state. 

Intriguingly, mixing of jet and envelope imprints on the 
jet a tight correlation (the "mixing attractor" of eq. 1461 
and figure IBJ between its Lorentz factor and energy per unit 
mass. This leads to a correlation between final opening angle 
and Lorentz factor, which in turn may be related to the lag- 
luminosity and variability-luminosity correlations observed 
in burst catalogs. Excessive mixing, however, has the effect 
of filtering rapid fiuctuations from jets. 

This work was stimulated by a visit to UC Santa Cruz 
and by interactions with Stan Woosley, Andrew MacFadyen, 
Alex Heger, and Weiqun Zhang while I was there. It was fur- 
ther motivated by conversations with and encouragement 
from Roger Blandford, Re'em Sari, Sterl Phinney and Josh 
Bloom during a visit to Caltech. It is a pleasure to thank 
Woosley, Blandford, and Phinney for their hospitality during 
those visits. I am especially grateful to Nicole Lloyd-Ronning 
for explaining the GRB luminosity-variability relation and 
suggesting clarifications. I also thank Andrei Beloborodov 
and Chris Thompson for discussing jet-envelope interactions 
and Chris Fryer for discussing progenitor scenarios. Com- 
ments from Beloborodov, Charles Dermer, Jonathan Tan, 
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Figure 7. Constraints on possible stellar progenitors for the case of GRB 000418 (assuming a ballistic jet and e-y = 10%). Stars in the 
black shaded region would collapse in the seven-second intrinsic duration of the burst. Those within the hatched region are excluded 
because they are too dense for a jet-cocoon structure to exist given the luminosity of this burst, and would develop a spherical blastwave 
instead. Those in the dotted region are also excluded because the GRB jet would take much longer than the observed duration to traverse 
their envelopes. 
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